Thermodynamics of strongly disordered spin ladders 
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We study antiferromagnetic two-leg spin- 1/2 ladders with strong bond randomness, using the real 
space renormalization group method. We find the low-temperature spin susceptibility of the system 
follows non-universal power laws, and the ground state spin-spin correlation is short-ranged. Our 
results suggest that there is no phase transition when the bond randomness increases from zero; for 
strong enough randomness the system is in a Griffith region with divergent spin susceptibility and 
short-range spin-spin correlation. 



I. INTRODUCTION 

One-dimensional (ID) spin systems are of interest ta 
physicists since the early days of quantum mechanics.EJ 
Considerable effort has been devoted to the theoretical 
study of antiferromagnetic (AF) spin chains, where some 
of the very few exact solutions _of interacting Hamilto- 
nians in physics were obtained,EH3 and remarkably rich 
low-energy physics we*e uncovered using various non- 
perturbative methods.cl Interest in these systems were 
also enhanced by the [-recent experimental realization 
of such model systemsp due to technological advances. 
More recently, considerable attention has focused on a: 
other class of ID spin systems, namely AF spin ladders.l 
These systems are made of two or more coupled spin 
chains. The physics of such systems are closely related 
to, but even richer than the spin chain systems. Fur- 
ther motivation for study of these systems comes from 
the similarity in structure between these systems and un- 
doped cuprates, and the discovery of superconductivity 
in them once charge carriers are introduced via doping. 

The ubiquitous randomness is known to have partic- 
ularly strong effects in low-dimensional systems. Re- 
cently, there has been rather extensive theoretical stud- 
ies of effects of disorder in spin chains. Most of these 
studies are based on the celebrated real space renor- 
malization group (RSRG) method introduced by Ma, 
Dasgupta and Hu (MDH) inpthe study of AF spin- 1/2 
chain with bond randomness ,cl and Bhatt and Lee in thg 
study of magnetic properties of doped semiconductors .0 
This method was elaborated and extended in great de- 
tail by Fisherja and applied (often with nontrivial ex- 
tensions) by a number of rf^pirStf^fWfhf 9 ] various dis- 
ordered spin chain modelsBE20E3^™ A variety of 
disorder-dominated phases have been found, whose low- 
energy physics is qualitatively different from their disor- 
der free counterparts. While the quantitative accuracy of 
the RSRG relies on the presence of strong randomness, 
it has been shownEl that even if the strength of random- 
ness is weak, it tends to grow as the RSRG proceeds to 
lower and lower energy scales, thus giving qualitatively 
correct (and sometimes asymptotically exact) low-energy 
behavior. Indeed, many predictions of RSRG have been 
confirmed by complementary analytical and numerical 
studies using other methods. 



Comparatively speaking, relatively few studies have fo- 
cused on effects of randomness on spin ladders. Effects 
of doping by non-magnetic impurities (or site .dilution) 
have been studied using quantum Monte CarloEJlia and 
mapping to Dirac fermions with random mass." 9 !! 20 ! The 
stability of the pure ladders against various types of weak 
randomness has been studied by Orignac and GiamachiEII 
using bosonization. In the present paper we study a two- 
leg AF spin- 1/2 ladder with strong bond randomness, 
using the RSRG. We believe our work is complementary 
to the previous studies, as the effects of bond random- 
ness and site dilution are quite different, and the RSRG 
is particularly suitable for studies of systems with strong 
randomness. 

While— .the present work was being completed, a 
preprinted appeared on the cond-mat archive, in which 
the authors used the RSRG as well as the density ma- 
trix renormalization group to study various disordered 
ladder models. While our work certainly overlaps with 
theirs, there exist two major differences, (i) Ref. ^ fo- 
cuses mainly on the distribution of the gap separating the 
ground and first excited states in finite clusters, while we 
study mainly thermodynamic properties and spin-spin 
correlation functions, (ii) Ref. |2^ has studied finite-size 
ladders with length up to 512. In our work we have stud- 
ied ladders with length up to 20,000, nearly a factor of 
40 bigger. The larger size is crucial to us for obtaining 
low-temperature, large distance behavior of the thermo- 
dynamic quantities and spin-spin correlation functions 
respectively. We will compare our results with those of 
Ref. ^ and previous studies whenever appropriate. 

Our main results are summarized as follows. We 
find the thermodynamics of the two-leg spin ladder re- 
mains non-universal, and the spin-spin correlation re- 
mains short-ranged, even in the strong-randomness limit. 
This is very different from the random— A F spin- 1/2 or 
spin-1 chain, wheperweak (for spin-l/2&E3) or sufficiently 
strong (for spin-lc3c3) randomness drive the system into 
the Random Singlet (RS) phase with universal thermo- 
dynamics and power-law spin-spin correlation. For suffi- 
ciently strong randomness, the spin susceptibility of the 
ladders exhibits power-law divergence as the tempera- 
ture T -> : x(T) - , where (3 varies continuously 
with randomness strength. Combined with short-range 
spin-spin correlation, this is characteristic of quantum 
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Griffith behavior. Such behavior appears to persist even 
when the interchain coupling of the ladder is made very 
weak, suggesting weak interchain coupling immediately 
destabilizes the RS phase that controls the low-energy, 
long-distance physics of the decoupled chains. 

The remainder of the paper is arranged as follows. In 
section II we introduce the model Hamiltonian we study, 
briefly review the RSRG method and its application to 
random spin chains, and discuss the necessary extensions 
we need to make in order to apply it to the ladder system. 
In section III we present our numerical results. In sec- 
tion IV we discuss the implications of our results, make 
contact with related theoretical and experimental work, 
and state our conclusions. 



II. MODEL AND THE 
RENORMALIZATION-GROUP SCHEME 

Consider an antiferromagnetic nearest-neighbor 
Heisenberg spin-I/2 ladder. The Hamiltonian for a 
two-leg ladder is given by: 



induce an effective coupling between spins 1 and 4 of the 
form: 



H 



N-l 

£ 

i=l, j=l,2 



N 

^i^id ' + KiSi^i ■ Si,2, (1) 



where N is the number of spins on a single chain, Sij is 
a spin-I/2 operator, and the positive coupling constants 
Jij (couplings along the chains, or legs of the ladder) and 
Ki (couplings between the chains, or along the rungs of 
the ladder) are distributed randomly according to some 
probability distributions P\\(Ji,j) and P±(Ki). 

In this work we use_the real space renormalization 
group (RSRG) methodQ to study the Hamiltonian (1). 
We begin by briefly reviewing its application to the ran- 
dom AF spin- 1/2 chains to illustrate the basic ideas be- 
hind it. In this approach one identifies the strongest cou- 
pling of the system, say, J2 that couples spins 2 and 3, 
and the two neighboring spins that are coupled to this 
strongly coupled pair. The Hamiltonian of this four-spin 
cluster is given by 



where 



H — H + Hi, 



Hq — J2S2 ■ S3, 

Hj = J1S1 ■ S2 + J3S3 ■ S4. 



(2) 



(3) 



In the presence of strong randomness, J2, being the 
strongest coupling in the system, is likely to be much 
stronger than other two couplings J\ and J3. In this case 
to a very good approximation spins 2 and 3 form a sin- 
glet pair in the low-energy states of the entire system and 
become inert. The weak perturbation Hi to this pair in- 
duces virtual transitions to the excited (or triplet) states 
of Hq; the main effect of such virtual transitions is to 



H, 



eff 



to second order in Hi, J is given by : 



J 



J1J3 
2J 2 



> 0. 



(4) 



(5) 



In essence, the RSRG procedure replaces the four spin 
cluster by spins 1 and 4, which are the active degrees of 
freedom at low energies, with a new effective AF bond J 
that couples them. J is typically much weaker than the 
original bonds (Ji, J2 and J3) so the bond distribution 
broadens and the energy scale is lowered. The decimation 
does not change the lattice structure of the chain, as after 
spins 2 and 3 are decimated, spins 1 and 4 become near- 
est neighbors, thus the new effective Hamiltonian still 
describes an AF spin-1/2 chain with nearest-neighbor 
interactions, and this procedure can be repeated. The 
renormalization scheme is depicted in Fig. [l]a. 

When we apply this technique to the ladder systems, 
some new features not encountered before appear, and 
corresponding modifications to the spin decimation pro- 
cedure described above need to be made. Firstly, the 
structure of the ladder gets distorted as soon as RSRG 
is applied, in contrast to the chain case. This situation 
requires us to keep track of the structure of the system 
as we decimate the spins and bonds. Secondly, when 
one decimates strongly coupled spin pairs, both AF and 
ferromagnetic (F) effective bonds are generated; these F 
bonds can lead to effective spins higher than spin-1/2 
at low energies. The initial renormalization step for the 
ladder is illustrated in Fig. [l]b, from which we can see 
how the ladder structure is distorted and ferromagnetic 
interactions are generated. These generated F bonds are 
much weaker than the original bonds that get decimated. 
However as we move on, more and more spins get deci- 
mated and the energy scale is lowered so at some point 
the generated ferromagnetic bonds can become impor- 
tant because it might be the strongest bond in the sys- 
tem. This spin pair, instead of forming a singlet, form a 
triplet or an effective spin-1 object at low-energy. Clearly 
larger spins can also be generated by RG as energy scale 
gets lowered. This situation is illustrated in Fig. Qc. In 
short, we need to keep track of both the lattice struc- 
ture, and size of the spins, together with the coupling 
constants (which can now be either AF or F) in our RG 
procedure. 

Now we turn the discussion to some technical details 
involved in the application of RG scheme to ladder. First 
consider two spins connected by a strong AF bond. These 
two spins are also coupled to the other 4 spins as shown 
in Fig. [j]b. We will make a slight change of notation 
for our discussion here just for simplicity. We label the 
spins participating in the process by number 1-6. The 
Hamiltonian for the 6-spin cluster is given by 



H = Hn 



Hr 



(6) 
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where 

Ho = J23S2 ■ S3, 

Hi = J12S1 ■ S2 + J34S3 ■ S4 + J25S2 ■ S5 + J36S3 ■ Sq, 

(7) 

where is the bond between Si and Sj. This 6-spin 
problem is quite complicated to solve but it can be sim- 
plified using the fact that we can treat Hj as a pertur- 
bation to Ho, especially when the randomness is strong. 
It is easy to see that to second-order, Hi only generates 
pairwise interaction among the spins. It is thus only nec- 
essary to include a pair of spins coupled to the two spins 
connected by the strongest bond, when we consider the 
effective interaction between them. This fact simplifies 
the calculation as we can now, reduce a 6-spin cluster 
problem to six 4-spin clusterscil which can be classified 
into three different types of 4-spin clusters as represented 
by: spin 1, 2, 3, and 4 (1234), (1235), and (2356). The 
Hamiltonians for these clusters are given by : 

Hi = J23S2 ■ S3 + J12S1 ■ S2 + J34S3 ■ £4, 
H2 = J23S2 ■ S3 + J12S1 ■ S2 + J25S2 ■ S$, 

H3 = J23S2 ■ S3 + J25S2 ■ S5 + J30S3 ■ Sq. (8) 

Hi has the same form as Eqs.(|J) and (^) which lead to 
the recursion relation Eq.(|). Cluster 1235, given by the 
Hamiltonian Hi, is a new cluster structure not encoun- 
tered in the chain case. Second order perturbation cal- 
culation shows that there is a new negative effective in- 
teraction between S\ and 5 5 given by : 



j _ J12J25 n 



(9) 



i.e. we have a ferromagnetic interaction. Physically this 
is due to the fact that Si and £5 are both coupled anti- 
ferromagnetically to S 2 ; this makes it favorable to have 
Si and S5 parallel to each other, thus an effective ferro- 
magnetic bond is generated. Cluster 2356, given by the 
Hamiltonian H3, looks almost the same as cluster 1234 
except that spin 5 5 and Sq are already connected by an 
original bond J56. This original bond will be renormal- 
ized when J 23 is decimated together with J25 and J36. 
The renormalized bond is then given by : 



J: 



50 



J; 



56 



J25J36 

2J 23 ' 



(10) 



The generated interaction between S5 and Sq is anti- 
ferromagnetic because they are sitting on the opposite 
sub-lattices. 

As discussed earlier, the effective F bonds generated 
by RSRG can lead to formations of effective spins with 
size larger than 1/2. We thus need to incorporate this 
possibility in our scheme, and generalize the Hamiltonian 
in Eq.(|^) and (^) by giving arbitrary sizes to the spin 
operators in the Hamiltonian and by having either sign 
for the couplings. We treat Hi as a perturbation to Hq 



as before. In the space of degenerate ground states of Ho, 
the spins S 1 and S2 form a state of maximum total spin 
5 = S2 + S3 for ferromagnetic ( J 2 < 0) or of minimum 
total S = |5 2 — 53 1 for antiferromagnetic (J 2 > 0) while 
the spins Si and S4 can point in any direction. The 
degenerate ground states span the Hilbert space TL which 
is the product space of the spin spaces Si, S, and S4. 
Hi will partially lift the degeneracy in TL and induce an 
effective Hamiltonian in TL. The effective Hamiltonian 
can be calculated using the projection theorerrcj: 



H 



eff 



PHP. 



(11) 



where P is the projection operator that projects the 
full Hamiltonian H into the subspace where S is max- 
imum(minimum). The detail of this calculation is avail- 
able in Ref. Here we just give the final result. After 
the strong bond is decimated, we can write down the 
effective Hamiltonian H e f f as : 



Heff = J1S1 ■ S + J3S ■ Si + constant, 



(12) 



where 



, S(S + 1) + 5 2 (5 2 + 1) - S 3 {S 3 + 1) j 



25(5 + 1) 
5(5 + l) + 5 3 (5 3 + l)-5 2 (5 2 + l) 
25(5 + 1) 



Ja, (13) 



where 5 = 5 2 ± S3I depending on the sign of J 2 . 

In the case where J 2 > and 5 2 = 53, the ground state 
of the strong bond is a singlet and there is no effective 
spin left after decimation. Second-order perturbation ex- 
pansion yields a non-zero interaction between 5i and 54: 



J = |5 2 (5 2 + l)^. 

o J 2 



(14) 



It can be shown that the cases discussed above exhaust 
all possible situations we may encounter when applying 
the RSRG to a spin ladder. 

In implementing the RSRG procedure outlined above, 
one finds that each spin is coupled to more other spins as 
more and more spins are decimated, and the couplings 
can be either F or AF. There is, however, one major 
simplification due to the bipartite nature of the original 
lattice, which is also of physical importance as we discuss 
below. In the beginning we have a lattice structure which 
can be divided into two sublattices (A and B) in which 
spins in sublattice A get coupled only to spins in sublat- 
tice B. As we run our RG procedure this is no longer true. 
Not only spins from the same sublattice can get coupled 
together but also the sizes of the spins are no longer the 
same. It becomes a relevant question to ask where to put 
an effective spins formed by two spins with different sizes 
and what the types of interactions are between this effec- 
tive spin with the rest of the lattice. We apply the major- 
ity rule in our RG scheme to incorporate this situation. 
The idea of this rule is to put the effective spin formed 
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1 2 3 4 J23 J12 J34 J14 



12 3 S e ff J23 J12 J12 



A A B A + 

A A B B + 

B A B A + 

B A B B + 



+ 



+ 



+ 



+ 



TABLE I: Some possible sublattice combinations for S2 = S3 
and J23 > 



by two spins with different sizes connected by AF/F cou- 
pling on the sublattice where the larger spin is. Using 
this rule we are able to show that two spins sitting on 
opposite sublattices will always have AF couplings while 
those sitting on the same sublattice will always have F 
couplings. This is clearly true in the beginning; we will 
show below the RG procedure combined with the major- 
ity rule preserves this structure. Physically this simply 
reflects the fact that the nearest neighbor AF couplings 
on a bipartite lattice has no frustration; they prefer the 
spins in the same sublattice to be parallel, and in oppo- 
site sublattices to be antiparallel. 

Let us elaborate this idea in more detail to better un- 
derstand the majority rule. We have seen in our dis- 
cussion above that there are three different cases which 
exhaust all the possible combinations encountered in our 
RG procedure. First, we have two spins with the same 
size connected by AF coupling. Second, two spins with 
the same or different sizes connected by F coupling and 
third, two spins with different sizes connected by AF cou- 
pling. These three cases are shown in Fig. [l]b and c. In 
the first case we do not have to worry about applying the 
majority rule because there is no effective spin formed. 
The configuration is shown on as cluster 1234 on Fig. 
[I|b. We just use the recursion relation derived in Eq. (|5|) 
to determine the type of interactions between the spins 
which were the third nearest neighbors. Some of the pos- 
sible sublattice combinations for this case is shown in Ta- 
ble (§) . Here it is clearly shown that two spins sitting on 
opposite sublattices will have AF interactions and those 
sitting on the same sublattice have F interactions. 

The configuration for the second case is shown on Fig. 
|l|c. We have already seen from Table (Q) that for two spins 
to have a ferromagnetic coupling, they must be sitting on 
the same sublattice. In this case there is no ambiguity 
where to put the effective spin. We can choose the effec- 
tive spin to be located on the site where either S2 or S3 
is used to be located. We can figure out the sign of the 
renormalized couplings in the same way as it is done in 
Table (0). The renormalized coupling is given by : 



Jn = 



So 



S3 



-J 



12- 



(15) 



With this recursion relation and majority rule, we can 
determine the sign of the renormalized couplings for all 
combinations possible. This is shown in Table (||). The 
conclusion that AF coupling is always on opposite sub- 
lattices and F coupling is always on the same sublattice 



ABB 
B B B 



B 
B 



TABLE II: Possible sublattice combinations for and J23 < 0. 
The column S e ff gives us the sublattice where we put the 
effective spin. 

12 3 Seff J23 J12 J12 



A A B 
BAB 



A 
A 



+ 
+ 



+ + 



TABLE III: Possible sublattice combinations for S2 > S3 and 
J23 < 0. S e f f gives us the sublattice where we put the effec- 
tive spin. 



remains valid. 

The last case is when S% ^ S3 and J 23 > 0. The 
majority rule tells us to put the effective spin on the 
sublattice of the spin with bigger size. If S2 > S3, we 
put the effective spin on the sublattice in which S2 is 
sitting and vice versa. The recursion relations for the 
couplings are given by the equation : 



J12 — Jn 



S 2 + l 
S2 — S3 + 



1 



(16) 



The type of interaction between the effective spin and 
the rest of the lattice is shown in Table ( [II ) , where we 
take an example 5*2 > S3. The result is the same as 
the two previous cases where AF coupling is always on 
opposite sublattices and F coupling is always on the same 
sublattice. Should we change S2 < S3, the result would 
remain valid. Table (III) shows the configurations where 
S 2 > S 3 . 

We have thus shown that the application of the ma- 
jority rule will preserve the type of interactions between 
spins sitting on opposite sublattices or the same sublat- 
tice. If the spins are sitting on the opposite sublattices, 
the interaction is always antiferromagnetic and if on the 
same sublattice, the interaction is always ferromagnetic. 
This conclusion can be generalized to higher dimensions 
as long as the original AF interactions couple only spins 
sitting on opposite sublattices. 



III. NUMERICAL RESULTS 

We have carried out the renormalization scheme for 
the ladder as described in section II numerically, with 
length of the ladders up to 20,000. In the decimation 
process, we pick up the strongest bond, as defined by the 
absolute value of the bond strength,c£l decimate it, and 
calculate the renormalized couplings to the neighboring 
spins. This procedure is iterated until the number of 
spins in the ladder is about 3% of the original number of 
spins. The initial distributions are taken to be in power- 
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law form : 

F||(J 4J ) = (1 n ;./, ;•.()• ./,._, ■ 1: 

= ^ K r^<K i <A. (17) 

Here < a < 1 is the measure of disorder (the bigger 
a, the stronger the randomness strength), and < A < 
1 is the anisotropy parameter; in the limit A — > the 
two chains decouple. We use a power-law form for our 
initial distributions because in the case of random spin 
chains, fixed point distributions at low energies typically 
have a power law form; we can thus hope to be able to 
approach the low-energy fixed points faster by starting 
with a power law distribution. 

As discussed earlier, due to the presence of F bonds 
generated by RSRG, effective spins with size bigger than 
1/2 appear at low energies. One might think that such 
larger spins may proliferate, and the typical size of the 
spins may grow indefinitely, leading to to a phase domi- 
nated by weakly coupled large spins. This was found to 
be the case in spin chains with random AF and F cou- 
plings studied by Westerberg et al£3 We find, however, 
this is not the case in the present problem. We address 
the issue of proliferation of F bonds and large spins in 
Fig. ||, where data for a = and a — 0.6 (both with 
A = 1) are shown. We plot the ratio of the numbers of 
AF bonds and F bonds as a function of bond strength 
cutoff fi in (a). At the early stages of RG the system 
consists of a large fraction of AF bonds and a small per- 
centage of F bonds generated by the decimation process. 
As the energy scale is lowered more F bonds are gener- 
ated and more AF bonds are removed so the ratio of the 
number goes down. In the low energy limit, we find the 
number of F bonds is very close to the number of AF 
bonds. This can be seen more clearly at the insets in 
Fig. ||a. Even though the numbers of AF and F bonds 
are almost equal, the strengths of AF and F bonds be- 
have completely differently in this limit. AF bonds al- 
ways dominate the system. In (b) we plot the ratio of 
the average strength of AF and F bonds. When the bond 
cutoff £1 goes below 0.2, where the numbers of AF and F 
bonds are almost equal, the ratio of the averages grows 
rapidly which means the AF bonds are much stronger 
than the F bonds in the low energy limit. In (c) we plot 
the difference of the averaged logarithms of AF and F 
bonds; the exponential of this quantity reflects the ratio 
between typical AF and F bonds. Similar to (b), here we 
see the the difference grows very fast, again showing the 
dominance of AF bonds over the F bonds. In (d) we plot 
the sample averaged ratio of the number of spins larger 
than 1/2 to the total number of spins. Here we see that 
while larger spins do appear, their percentage remains 
small, and the percentage decreases with the cutoff f2 
going down in the low-energy limit. Another piece of in- 
formation that is not included in the figure is that most 
of the larger spins are spin-l's, with a very small per- 
centage of spin 3/2 and spin 2. We have not found any 
trace of spins larger than 2 in our simulations. We find 



qualitatively similar behavior in all initial distributions 
we have looked at, indicating this is generic. 

Physically, such behavior has its origin in the bipartite 
nature of the lattice structure of the 2-leg ladder. As we 
have shown earlier, the effective couplings generated by 
the RSRG is always AF between spins of opposite sublat- 
tices, and F between spins of the same sublattice. Since 
the numbers of spins in the two sublattices are the same, 
the number of F and AF bonds become very close in the 
low-energy limit. On the other hand spins in opposite 
sublattices tend to be closer to each other, leading to the 
fact that AF bonds dominate F bonds in strength. This 
in turn suppresses formation of large spins. 

Our most important results are presented in Figs. 3-6, 
where we plot the temperature dependence of the spin 
susceptibility, and the ground state spin-spin correlation 
function. The susceptibility is calculated as the following. 
We proceed with the RSRG until the bond cutoff fl is 
equal to the temperature T. We neglect contributions of 
spins that have already been decimated, and treat the 
remaining spins as free spins, thus their contribution to 
the susceptibility is just the Curie susceptibility. This is 
a good approximation as long as the bond distribution is 
broad. The total susceptibility is thus given by : 

xtot = ^J2 N ^ s + 1 ^ ( 18 ) 

^ s 

where N s is the number of spins left at energy scale Q = T 
for a given spin size s and the summation runs over all 
possible spin sizes. In Fig. |^ we plot the susceptibility per 
spin for different samples as a function of temperature for 
different disorder strength a, all with isotropic coupling 
(A = 1). In all cases we find the low-T susceptibilities can 
be fit quite well to power-law dependence on T: x ~ T - ' 3 ; 
the power-law exponent (5, which we obtain from a least- 
square fit to the low-T part of the data, is non-universal; 
it can describe both divergent (/? > 0) x f° r stronger 
randomness (larger a), or vanishing x {P < 0) for weaker 
randomness (smaller a), as T — > 0. It is worth noting 
that for A = 1, we always have (3 < a, and such behavior 
persists for very strong disorder like a — 0.9. Such be- 
havior is very different from random AF spin- 1/2 chain 
with any amount of randomness, or random AF spin- 
1 chain with sufficiently strong randomness, where the 
system flows to the so-called random singlet (RS) fixed 
point, in which the bond distribution is infinitely broad, 
the spin-spin correlation follows a universal power-law, 
and the susceptibility diverges in a universal manneiu: 

X ~l/(Tln 2 (f!o/T)). (19) 

Instead, the fact that we find power-law exponent to be 
always less than 1 indicate the width of the bond distri- 
bution is finite. Of course, in principle we can not com- 
pletely rule out the possibility that our system size (and 
correspondingly, temperature range) is not wide enough 
for us to approach the true low-T asymptotic behavior 
of x, which for strong enough randomness may be con- 
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trolled by a fixed point with infinitely broad bond distri- 
bution and universal. We believe, however, this is highly 
unlikely for the following reasons, (i) Our power-law fit 
already extends to a very wide range in T. In partic- 
ular, for a = 0.9, a single power-law fits all the data 
very well that is over eleven orders of magnitudes in T, 
with no indication of crossover to other behavior at low 
T. (ii) As we will see later, the spin-spin correlation 
function appears to be short-ranged, indicating that the 
long-distance, low-temperature physics is not controlled 
by a single scale- invariant fixed point. 

In the absence of interchain coupling, the ladder be- 
comes two decoupled random AF spin- 1/2 chains, where 
the long-distance, low-temperature physics is controlled 
by the RS fixed point and universal. To address how the 
system crosses over from one behavior to another we have 
studied how the susceptibility varies with the anisotropy 
parameter A. In Fig. ^ the susceptibility per spin for 
different values of A is presented, for a fixed a = 0.6. 
Again, we find non-universal behavior here. As we vary 
A from I to 0, the power-law exponent of the suscepti- 
bility increases continuously. In the case of A = we 
have decoupled chains and the susceptibility is expected 
to follow Eq. (19). While for a finite range of T it can 



be fit reasonably well to a power-law with f3 very close 
to 1, the small upward curvature of the data indicates (3 
would increase as one goes to lower T, consistent with 
Eq. ([H]). On the other hand a very weak interchain 
coupling (e.g., A = 0.025) leads to a significant change 
it j3, and there is no longer obvious upward curvature in 
the data. This suggests that a weak interchain coupling 
immediately destabilizes the RS fixed point. 

We now turn the discussion to the ground state spin- 
spin correlation function along the chain: 



• = (-!)*-' «S itk -S j<k », 



(20) 



where <<>> stands for both quantum and disorder av- 
erage. We calculate g(\i — j\) in the following way. We 
run the RSRG until all spins are decimated, and then 
simply count the number of singlet pairs formed for a 
given distance \i — j\, divide this number by the total 
number of pairs and multiply -the result by 3/4. In the 
RS phase, g(\i —j\) ~ \i — j|~ 2 .tl In Fig. [|a we study how 
the interchain interaction affects the correlation along the 
chain by varying the anisotropy parameter A. Fitting the 
data to a power-law dependence: g(\i — j\) ~ \i — , 
we obtain v = 1.97 for A = (decoupled chain ease), 
which is very close to the analytical result v = 2.0 For 
nonzero A, the correlation decays much faster than that 
of the chain. Even a small amount of interchain inter- 
actions (say, A = 0.001) change the behavior of the cor- 
relation considerably. We can see a downward curvature 
in the data, which is particularly obvious for A = 1 and 
0.5, indicating the short range (decaying faster than any 
power-law) behavior of the correlation. If we try to fit 
the ground state correlation for non-zero A to a power- 
law, we would get considerably larger power-law expo- 
nent is, even for A as small as 0.001. This strongly sug- 



gests that introduction of interchain interactions imme- 
diately destabilizes the RS phase that controls the low- 
energy of the decoupled chains, and leads to short-range 
spin-spin correlation in the ground state. We have also 
calculated how the correlation changes as we vary a for 
fixed A = 1, in Fig. ^3. Here we find while stronger 
randomness (larger a) tends to enhance correlation at 
large distances, the correlation is still short-ranged for 
very strong randomness (a = 0.9) as evidenced by the 
downward curvature of the data. 

One general concern in numerical calculations of the 
kind discussed here is finite-size effects. We show in Fig. 
^ that the system sizes we use in this work are large 
enough that the finite-size effects are negligible. The 
sample averages of the susceptibility per spin do not show 
any noticeable fluctuations as the system sizes are varied 
from N = 2, 000 to N = 20, 000. The same is true for 
the ground state spin-spin correlation. At large separa- 
tion there are some variations due to sample to sample 
fluctuations. We can thus safely conclude that the finite- 
size effect is negligible in our study. 



IV. SUMMARY AND DISCUSSIONS 

In this work we have used the RSRG method to study 
an AF two- leg spin- 1/2 ladder, with strong bond ran- 
domness. We find that the spin susceptibility is non- 
universal, and the ground state spin-spin correlation is 
short-ranged, for any randomness and interchain cou- 
pling strength. For sufficiently strong randomness or suf- 
ficiently weak interchain coupling, the spin susceptibility 
exhibit power law divergence as T — > 0, which is charac- 
teristic of quantum Griffith behavior. 

Melin et alr3 used the RSRG method as well as den- 
sity matrix RG to study the distribution of the energy 
gap separating the ground and first excited states in clus- 
ters (with length up to 512) of AF two-leg spin-1/2 lad- 
ders. They find that the dynamic exponent z that char- 
acterizes this distribution is non universal and depends 
continuously on randomness strength. Based on this they 
conclude that the low-energy physics of the system is con- 
trolled by a fixed point with a finite width in the bond 
distribution function, and the system is in a quantum 
Griffith phase. Our results and conclusions agree with 
theirs. 

It is by now well established that in the absence ran- 
domness, the two-leg AF spin-1/2 ladder supports a finite 
excitation gap, and the spin-spin correlation is short- 
ranged. It is generally true that randomness tends to 
introduce low-energy excitations, which can lead to di- 
vergent spin susceptibilities as found here. Our results 
indicate however, despite the low-energy excitations in- 
troduced, the phase with short range spin-spin correla- 
tion appears to be stable against any amount of random- 
ness. This is certainly consistent with Ref. |2l], where 
the authors find the pure ladder to be remarkably stable 
against various kinds of disorder. On the other hand this 
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is very different from the AF spin-1 chain, where suffi- 
ciently strong bond randomness drives the system from 
the Haldane phase to the random singlet phase with uni- 
versal thermodynamics and power-law spin-sa in , c prrela- 
tion, through a second-order phase transition.li3liJ 

As discussed earlier, for the present system the bond 
distribution has a finite width in the low-energy limit, 
no matter how strong the randomness is initially. This 
indicates that the RSRG method is not asymptotically 
exact when applied to the present model. However this 
method is quantitatively accurate as long as the random- 



ness is strong, we thus believe the qualitative conclusions 
we draw from our results are robust. 
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(c) 



FIG. 1: (a) The renormalization scheme for a four-spin problem with AF couplings, as encountered in random AF spin-1/2 
chains. Here the strongest bond J2 is decimated, together with the neighboring bonds Ji and J3, yielding an effective interaction 
J between what were the third-nearest neighbors, (b) Schematic diagram for decimation in ladder. The dashed lines are the 
renormalized couplings. The thick dashed lines are the ferromagnetic couplings generated in the decimation process, (c) The 
renormalization scheme for four spin problem where the strongest bond is ferromagnetic. The two spins connected by F coupling 
forms an effective spin object having renormalized interactions with its neighbors. 
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FIG. 2: Numerical simulation results of the proliferation of ferromagnetic bonds and larger spins, for two different initial bond 
distribution (A = 1 in both cases), (a) The ratio of numbers of AF and F bonds; (b) the ratio of average strength of AF and 
F bonds; (c) the difference of the averaged logarithm of AF and F bond strengths, and (d) the ratio of number of spins larger 
than 1/2 and the total spins, all as functions of energy scale fi. All samples used have size N = 20, 000 and in (a)-(c) different 
symbols represent different samples. 
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(°) kg T (d) kg T 

FIG. 3: The susceptibility per spin as a function of temperature for several disorder strengths a. The anisotropy parameter 
A = 1 is fixed, and the system size is N = 20, 000. The power law exponents, (3, are calculated using the least square fit to the 
low-temperature data. The different symbols in the figure correspond to different samples. We do not take the sample average 
because the sample to sample variations are small. 



11 






-in 4 
10 




10 J 




10 2 


""eh. 










10 1 








10° 




10" 1 



. 1 1 II 


, 1 1 '1 1 1 '1 1 ' '■ 

A = 0.025 




p = -0.712 - 






. , .1 





10" 5 10" 4 10" 3 10" 2 10" 1 10° 



(b) 



k B T 



10 J 



10* 



10' 



10 u 



10" 



1 — ' ,_ T — ' ,_ T — ' 

A = 0.05 
P = -0.670 




10"° 10" 4 10" 3 10" 2 10" 1 10° 
(c) k B T 



O) 



10 J 



10* 



10' 



10 u 



10" 



r ~ r l — ' ,_ T — ' ,_ T — ' ,_ T — ' 

A = 0.1 

^ % P = -0.621 



(d) 




jj , ■_j , , . .i 



k B T 



10" 5 10" 4 10" 3 10" 2 10" 1 10° 



CD 



10" F 



10' r 



10 u r 



10" 



10"' 



(e) 



"i — i — i — a 

A = 0.5 



10"' 



10" 
k B T 



P = -0.505 




10" 



10 u 




FIG. 4: The susceptibility per spin as a function of temperature for a given a — 0.6 with varying anisotropy parameter A. 
System size fixed to be N = 20, 000. The power law exponent increases as we decrease the anisotropy parameter. For A = 
we have decoupled chains. We plot all different samples in the figure without taking the sample average because the variations 
are very small. 
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FIG. 5: The sample-averaged spin-spin correlation along the chain (a) with varying A and a = 0.6, (b) with varying a and 
A = 1, for N = 20,000. The fit to power-law behavior yields the following power-law exponents : (a) v = 1.97 for A = 0, 
v = 2.56 for A = 0.001, v = 2.87 for A = 0.5, v = 3.23 for A = 1, and (b) v = 3.45 for a = 0.3, v = 3.23 for a = 0.6, v = 2.89 
for A = 0.9. 
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FIG. 6: (a) The susceptibility per spin as a function of temperature, (b) the spin-spin correlation along the chain, for a — 0.6 
and A = 1 with varying system sizes. No significant variation in these two quantities as the system sizes are varied. Both (a) 
and (b) are sample averages. 



